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Abstract. We predict mass-loss rates for the late evolutionary phases of low-mass stars, with special emphasis 
on the consequences for the morphology of the Horizontal Branch (HB). We show that the computed rates, 
as predicted by the most plausible mechanism of radiation pressure on spectral lines, are too low to produce 
EHB/sdB stars. This invalidates the scenario recently outlined by Yong et al. (2000) to create these objects by 
mass loss on the HB. We argue, however, that mass loss plays a role in the distribution of rotational velocities of 
hot HB stars, and may - together with the enhancement of heavy element abundances due to radiative levitation 
- provide an explanation for the so-called "low gravity" problem. The mass loss recipe derived for hot HB (and 
extreme HB, sdB, sdOB) stars may also be applied to post-HB (AGB-manque, UV-bright) stars over a range in 
effective temperatures between 12 500 - 40 000 K. 
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1. Introduction 

Over the last decades, both observational and theoretical 
efforts have been devoted to the investigation of the ob- 
served distribution of stars along the Horizontal Branch 
(HB) of galactic Globular Clusters (GCs). Although 
canonical stellar evolution theory has provided a gen- 
eral consensus on the evolutionary phase corresponding 
to the HB sequence, and convincingly demonstrated that 
its morphology is most strongly affected by cluster metal- 
licity (the 'first' parameter; Sandage & Wallerstein 1960), 
many problems remain. The most striking controversy in- 
volves the wide variety of HB morphologies among clusters 
with similar metallicities (the 'second parameter' problem; 
Sandage & Wildey 1967, van den Bergh 1967). Candidate 
second parameters are cluster age (e.g. Lee et al. 1994 and 
references therein) , mass loss along the Red Giant Branch 
(RGB) (Catelan et al. 2001 and references therein), ro- 
tation and deep helium mixing (Sweigart 1997), dynami- 
cal interactions involving binaries and even planets (Soker 
1998), as well as environmental effects in high-density en- 
vironments (Fusi Pecci et al. 1993). 

The identification of the second parameter is espe- 
cially relevant to the formation of Extremely blue HB 
(EHB) stars, which are thought to be responsible for 



the ultraviolet upturn phenomenon in elliptical galaxies 
(Greggio & Renzini 1990; Dorman et al. 1995). The pres- 
ence of EHB stars as blue 'tails' in clusters (Ferraro et 
al. 1998; Piotto et al. 1999), as well as sdB/sdO stars 
in the field (Greenstein 1971; Green et al. 1986), has 
inspired modern-day research to explain their formation 
both through mechanisms that produce high mass loss 
along the RGB (Soker et al. 2001 and references therein), 
as well as through binarity (Mengel et al. 1976, Heber et 
al. 2002). 

Further puzzles in HB morphology concern the issues 
of HB 'gaps' (Newell 1973) - specific regions along the 
branch that are significantly underpopulated^, and a rel- 
atively new, unexplained, but ubiquitous feature is the 
so-called Strbmgren u-jump at an effective temperature of 
T eff ~ 11 000 K (Grundahl et al. 1999), possibly coinciding 
with a jump in log g (Moehler et al. 1995), and an unex- 
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1 It has been claimed that the positions of the gaps along 
the HB in different galactic GCs are the same within current 
empirical uncertainties (Ferraro et al. 1998), however it is not 
clear whether these gaps mark regions with specific effective 
temperatures, or whether they correspond to constant mass 
loci (Piotto et al. 1999). Note that Catelan et al. (1998) have 
challenged the existence of gaps at the same positions in all 
clusters. 
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plained absence of fast rotators above this temperature 
(Behr et al. 2000; Recio-Blanco et al. 2002). 

As discussed by Grundahl et al. (1999), the Stromgren 
u-jump may be due to atmospheric diffusion by radia- 
tive levitation of heavy elements, as both Glaspey et al. 
(1989) and Behr et al. (1999) found striking abundance 
anomalies in blue HB stars, with iron enhancements of 
up to three times the solar value. Moehler et al. (2000) 
have shown that the enhancement of heavy elements in 
spectroscopic analyses may partially solve the problem of 
the anomalously low gravities along the blue HB, but the 
discrepancy is still present at the level of Alogg w 0.1 
dex for stars in the range 15 000 < T cS < 20 000 K. Even 
more so, the first two mentioned HB features - the blue 
tails and the gaps - are still an enigma^], and it is not 
at all obvious whether they originate from a mechanism 
working in a prior evolutionary phase (on the RGB) or if 
they are due to a process working 'in situ' once the star 
has settled on the HB. One of the options that may help 
in explaining the above-mentioned problems is mass loss 
on the HB. 

It is worth mentioning that the colour width of the 
hottest gap is so small that changes of the order of a few 
times 10 -3 Mq in the total mass are capable to move the 
star away from its initial location, far enough as to produce 
an underpopulated region in the H-R diagram (HRD). For 
this to occur one needs to identify a mass-loss mechanism 
which efficiency rapidly increases at the specific effective 
temperatures of the gap. 

The hypothesis that a mass-loss mechanism may be 
at work during the HB evolution was first entertained by 
Wilson & Bowen (1984). They suggested that an increased 
mass-loss efficiency, when crossing the RR Lyrae instabil- 
ity strip, could provide an explanation for the HB mass 
distribution in a more natural way than the alternative 
of a stochastic variation in the amount of mass lost dur- 
ing the prior RGB phase. The topic of mass loss on the 
HB was further addressed by Koopmann et al. (1994), but 
they concluded that constant mass loss in the RR Lyrae 
strip was incapable of providing an explanation for the 
HB mass dispersion or the RR Lyrae period change dis- 
tribution. Additionally, mass loss during the central He- 
burning phase was suggested by Michaud et al. (1985) and 
Bergeron et al. (1988) in order to explain the large silicon 
underabundances in some HB stars. More recently, Yong 
et al. (2000) performed accurate evolutionary computa- 
tions with mass loss, and suggested that mass-loss rates 
of the order of 10" 9 - 10" 10 Moyr" 1 for HB stars in the 
metal-rich cluster NGC 6791 can force these stars to move 
to a bluer position and thus lead to the production of EHB 
stars. If correct, this scenario could provide an explana- 
tion for the presence of extended blue tails along the HB 
of some metal-rich GCs, such as NGC 6441 and NGC 6388 
(Rich et al. 1997), although it would not be able to explain 

2 Note that Brown et al. (2001) and Sweigart et al. (2002) 
have provided a theoretical framework which could explain the 
hot gap in the HB of NGC 2808. 



the upward sloping of the HB in these clusters (Raimondo 
et al. 2002 and references therein). The main problem with 
the proposed scenario, however, is that no physical mech- 
anism for mass loss was proposed and that the adopted 
mass-loss rates were completely 'ad hoc', as there are nei- 
ther observational data indicative of mass loss on the HB 
available, nor any predictions. 

Our aim in the present paper is to alleviate current 
shortcomings by computing radiation-driven wind mod- 
els and mass-loss rates for low-mass blue stars, and to 
subsequently investigate their influence on HB evolution- 
ary models. Blue HB stars are located in a region of the 
HRD, where the stars are hot (with T e g between 10000 
and 35 000 K), and relatively bright, and radiation pres- 
sure forces can therefore be considered a natural driv- 
ing mechanism. Although there may be other processes 
that could possibly drive a wind, such as pulsations all 
other wind-driving options are much less well-understood 
than radiation pressure on spectral lines. Radiation-driven 
wind models have been developed in the 1970s by Lucy & 
Solomon (1970) and Castor et al. (1975). In more recent 
days, the models have been very successful in predicting 
the values observed in O supergiants (Vink et al. 2000). 
The direct application of these predictions to HB stars, 
such as the use of the mass-loss recipe provided by Vink 
et al. (2000) would however involve a rather large and 
dangerous extrapolation by four orders of magnitude in 
stellar luminosity. 

As far as the 'gaps' along the HB are concerned, 
radiation-driven wind models for OB supergiants predict 
that the efficiency of mass loss jumps strongly by a factor 
of five at spectral type Bl (Vink et al. 1999, 2000). This 
is close to the position where the evidence for a gap in HB 
morphology is strongest. A mass- loss rate of the order of 
10 -10 — 10~ n M Q yr _1 could be sufficient to explain the 
presence of the gap located at T e g~20 000 K; given an 
HB evolutionary timescale of w 10 8 years with mass loss 
at this rate leads to a total amount of a few times 10 -3 
Mq sufficient to move an HB star by wlOOO K, and so 
creating a 'gap'. 

The above-mentioned issues, i.e. the presence of EHB 
stars, gaps, and anomalous abundances in HBs and sdB 
stars, prompted us to compute radiation-driven wind 
models for HB stars; to predict mass-loss rates for these 
objects, and subsequently explore their influence on evolu- 
tionary models. The mass loss computations may also pro- 
vide valuable ingredients for HB angular momentum evo- 
lution and chemical separation calculations of sdB stars 
(see Unglaub & Bues 2001). 

The outline of the paper is as follows. In the next sec- 
tion we describe the approach used for computing mass- 

3 Recently a new class of variable stars has been discovered 
in the field, the so-called EC 14026 (Kilkenny et al. 1997), which 
have been identified as hot HB stars and their progeny. Even 
though no clear identification of similar variables in GCs have 
been obtained, there is no a priori reason for the lack of this 
kind of pulsation among cluster HB stars. 
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loss rates, as well as the assumptions adopted in the nu- 
merical computations; in Sect. 3 we discuss the results con- 
cerning the mass-loss efficiency, where the dependence of 
M on the main evolutionary parameters, the luminosity, 
effective temperature and stellar mass, as well as stellar 
metallicity, is presented. In Sect. 4, we provide an analyti- 
cal relation for M as a function of the quoted parameters, 
which is useful for computing the mass-loss rates in evo- 
lutionary computations, and we investigate the effects of 
our recipe on HB stellar evolution (Sect. 5). In Sect. 6, we 
study the implications of mass loss regarding the "zoo" of 
problems in HB morphology that occur for effective tem- 
peratures larger than ~ 10 000 K, in particular the effects 
of mass loss on rotational velocities and the logg jump. 
Final remarks and conclusions will close the paper. 

2. Method and Assumptions in computing M 

We compute mass-loss rates for HB stars under the hy- 
pothesis that radiation pressure on spectral lines drives a 
stellar wind on the HB. Although this does not imply that 
there are no other physical mechanisms operating during 
the HB phase that could possibly drive a wind, it is the 
most sophisticated wind theory known, and it has been 
very successful in explaining the observed mass-loss rates 
of hot massive stars. 

2.1. The Monte Carlo method to calculate M 

The description of the radiative wind driving with our 
method is based on a Monte Carlo technique that was 
first introduced by Abbott & Lucy (1985). This approach 
naturally accounts for photon-interactions with different 
metal ions, as the photons try to escape from the stellar 
wind. In the Monte Carlo model used here (mc-WIND, 
de Koter et al. 1997, Vink et al. 1999), the momentum 
deposition is calculated using the Sobolev approximation 
by following the fate of a large number of photons that are 
released from below the stellar photosphere. To obtain a 
consistent solution, several wind models are calculated to 
find the mass-loss rate that is consistent with the radiative 
acceleration (see also Lucy & Abbott 1993). 

The calculation of radiation pressure with this method 
requires the input of a model atmosphere. The model at- 
mospheres used in this study are the non-LTE unified 
Improved Sobolev Approximation code (isa-wind), which 
treats the photosphere and wind in a unified manner (dis- 
tinct from the so-called "core-halo" approaches). For de- 
tails of the code we refer the reader to de Koter et al. 
(1993,1997). The chemical species that are explicitly cal- 
culated in non-LTE are H, He, C, N, O, and Si. The iron- 
group elements are treated in a generalised version of the 
"modified nebular approximation" (Lucy 1987, 1999). 

2.2. Assumptions in the models 

The model depends upon the assumption that the plasma 
behaves as a single fluid. As long as a large number of 
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collisions between the accelerating (C,N,0, and Fe-group) 
and non-accelerating (H and He) particles ensures a strong 
coupling, one can safely treat the wind as a single fluid. A 
simple condition for this so-called 'Coulomb coupling' is 
given by Lamers & Cassinelli (1999, p. 193): 

< 5.9 x 10 16 (1) 

M 

where L„ is in Lq, v is in kms -1 and M is in M Q yr _1 . 
Although this Coulomb condition is easily satisfied for the 
dense winds of massive O stars, it may break down for 
weaker winds, such as the winds from main-sequence A 
and B stars (Springmann & Pauldrach 1992; Porter & 
Drew 1995; Babel 1996; Krticka & Kubat 2000). However, 
using typical values for an EHB star, logL* = 1.3; v — 600 
km s _1 ; and logM = —11.70 (see computations in Sect. ||), 
we find that the ratio (L*v/M) from Eq. (|l|) is ten times 
smaller than the quoted value of 5.9 x 10 16 , indicating that 
the condition of Coulomb Coupling is fulfilled. It remains 
yet to be seen if the use of the Sobolev approximation is 
valid for weaker winds (see Owocki & Puis 1999). 

2.3. The adopted parameters 

Effective Temperatures The models have effective temper- 
atures between 12 500 and 35 000 K with stepsizes of 2 500 
K. We have checked this choice of gridding by comput- 
ing additional models at intermediate temperatures, which 
showed our initial resolution to be entirely adequate. 

Luminosities and Masses The values for the stellar lumi- 
nosity (L) and mass (M) were taken from evolutionary 
models (Cassisi & Salaris 1997; Zoccali et al. 2000) for HB 
stars computed under different assumptions regarding the 
initial chemical composition. 

Metal abundances The adopted heavy elements distri- 
bution corresponds to the solar scaled one provided by 
Anders & Grevesse (1989), adopting Z Q = 0.019. For 
non-solar metallicity Z , the helium (Y) and hydrogen (X) 
abundances are adjusted in the following way: 

where the primordial helium abundance Y" p = 0.24 
(Audouze 1987) and (AF/AZ) = 3 (Pagel et al. 1992). X 
is then simply given by: 

X = 1 — Y — Z (3) 

Velocity Law We calculated M for wind models with a 
/3-type velocity law for the accelerating part of the wind: 

v(r) = Voo fl - -^J (4) 
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Fig. 1. Mass loss predictions for HB stars as a function of 
effective temperature. The solid line represents the com- 
putations for a solar metallicity. The dotted line is for a 
metallicity three times solar. The models are calculated 
for a constant mass of M = 0.5M Q and log(L/L Q ) = 1.5. 



Below the sonic point, a smooth transition from this ve- 
locity structure is made to the velocity that follows from 
the photospheric density structure. A value of (3 — 1 was 
adopted in the accelerating part of the wind. This is a 
typical value for OB supergiants (Groenewegen & Lamers 
1989, Puis et al. 1996), but future mass loss observations 
of HB stars should be able to show whether or not this is 
also an adequate description for these low mass objects. 
Note that it has been demonstrated that over a (3 range 
0.7 - 1.5 the M predictions are insensitive to the adopted 
value of (3 (Vink et al. 2000). We further assume a ratio 
«oo csc = 1-0, a s roughly indicated by Voo observations 
for a handful of SdO stars by Hamann et al. (1981) and 
Howarth (1987). 



3. The mass-loss predictions 



Using the procedure described in Sect. 2.1, we have calcu- 
lated mass-loss rates as a function of T e g with tempera- 
tures in the range between 12 500 and 35 000 K. This was 
performed for luminosities in the range log (L„/Lq) = 1.3 
-1.7 and masses in the range M* = 0.5 - 0.7 M & . 

3.1. Mass loss as a function ofT c g 

The results of our predictions of HB mass loss as a func- 
tion of effective temperature are presented in Fig. |l|. The 
solid line represents the computations for solar metallic- 
ity. To check the generality of this behaviour we have also 
computed mass loss as a function of T e g for somewhat dif- 
ferent input parameters, represented by the dotted line, for 
a metallicity three times solar. In both cases there appears 
to be a slight decrease of mass loss as a function of decreas- 
ing effective temperature. This can be attributed to the 
gradual shift of the flux maximum towards longer wave- 



Fig. 2. Mass loss predictions for HB stars as a function 
of stellar luminosity. These models are calculated for a 
constant mass of M = O.5M , r o ff = 20 000 K, and solar 
metallicity. 



lengths, and as the number of lines present in the spec- 
trum is smaller at higher wavelengths, the line acceleration 
decreases, reducing the mass-loss rate. Superimposed on 
this, one may have expected to see jumps, where M could 
increase due to recombinations of important line-driving 
ions. As mentioned earlier, for OB supergiants the mass 
loss increases steeply by a factor of five due to the recom- 
bination of Feiv to Fein at spectral type Bl. 

Figure || however, indicates that these so-called "bi- 
stability" jumps are absent for HB stars. This is probably 
due to the lower wind densities in HB stars in compar- 
ison to OB supergiants. In fact, these HB computations 
(for solar metallicity) are more comparable to the OB su- 
pergiant calculations at very low metallicities. Vink et al. 
(2001) have shown that at lower wind densities, the winds 
are no longer driven by Fe, but that the line driving by 
CNO-like elements takes over, and the dramatic recombi- 
nations are much less pronounced, or even absent at the 
temperatures under consideration here. 

We conclude, that bi-stability jumps are not present 
for HB stars and we apply a fit through our computed dat- 
apoints over the complete temperature range. This tem- 
perature dependence of mass loss is later incorporated into 
our mass loss recipe for HB stars (Sect.|j). 

3.2. Mass loss as a function of L, M and Z 

Predictions of HB mass loss as a function of other stel- 
lar parameters, namely L, M , and Z have also been per- 
formed. As an example, the results for mass loss as a 
function of L are shown in Fig.|^. The results of M as 
a function of L, M and Z do not yield any surprises. 
As expected, mass loss increases for increasing L and Z, 
but decreasing M. Although the behaviour is qualitatively 
similar to the O star recipe in Vink et al. (2000), the de- 
pendencies are somewhat different. The values of these de- 
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pendencies and mass-loss rates are discussed in the next 
section. 

4. Mass loss recipe 

In this section, we present a mass loss recipe for HB stars 
as a function of basic stellar parameters. To obtain the 
recipe we have determined four separate dependencies, 
and checked if they were independently applicable. As this 
was found to be the case, we have com bine d th e four in- 
dependent parameters from Sects. ^1 and ^2 and con- 
structed the following analytical relationship for HB mass 
loss: 



log M = - 11.70 (±0.08) 

+ 1.07 (±0.32) log(T off /20000) 
+ 2.13 (±0.09) (log£* - 1.5) 
- 1.09 (±0.05) log(M*/0.5) 
+ 0.97 (±0.04) log(Z„) 

derived for : 

12 500 < T cff < 35 000 K 
1.3 < logL, < 1.7 
0.5 < M* < 0.7 
0.1 < Z* < 10 



(5) 



where T c g is in Kelvin and L„, Z*, and M* are all given 
in solar units. Note that the computatio ns h ave been per- 
formed for the case Voo — v csc (see Sect. 2J5). As discussed 
in Vink & de Koter (2002), M is maximal for this value, 
and Eq. (||) therefore provides an upper limit to the mass- 
loss rate. In case v^ is not equal to w csc , one can use 
the ( Voo/ v C sc) dependence of M as derived in Vink et al. 
(2000) for OB stars. The root-mean-square (rms) differ- 
ence between the results following from Eq. (|J) and the 
actual model computations is 0.08 dex in log M. This im- 
plies that the HB formula yields a good representation of 
the actual HB mass loss computations. 

Alternatively, if we put the HB stellar parameters into 
the recipe of Vink et al. (2000) for massive O stars (and 
keep Voo fixed to v CS c)j we find an rms difference of 0.44 
dex in the mass loss, leading to the conclusions that an 
extrapolation of the O star recipe by four orders of magni- 
tude in stellar luminosity would have resulted in mass-loss 
rates that are systematically too high by about a factor 
of two. One may still wonder whether the HB mass loss 
recipe of Eq. (||) is also applicable to stars with stellar pa- 
rameters for which the recipe was not specifically derived. 
To check whether Eq. (||) may safely be used over a wider 
range in stellar parameters, we have performed mass loss 
calculations for the winds of post-HB stars and compared 
these actual calculations with the results from the recipe. 
We find that the mass loss recipe may be used for other 
classes of low-mass blue objects as well, as long as the 



desired accuracy is within a factor of ~ 2. As there are 
hardly any mass-loss predictions available for these types 
of objects either, the HB formula may be applied to all 
hot, low-mass stars, of the types: HB, EHB, sdB, sdOB, 
post-HB, AGB-manque, UV-bright stars, and extreme he- 
lium stars, as long as their effective temperatures are not 
significantly higher than ~ 40 000 K. This because exist- 
ing mass loss calculations, such as the ones presented here, 
but also those by Pauldrach et al. (1988) for Central Stars 
of Planetary Nebulae have the problem that line lists be- 
come incomplete with respect to higher ionisation stages. 
Note that we do not expect problems with extrapolating 
Eq. (H) to effective temperatures as low as « 8 000K. A 
computer routine of the HB mass loss recipe is available 
upon request or on the Web.Q 

5. Mass loss effects on HB evolution 

5.1. Assumptions in Evolutionary Tracks and Outer 
Boundary conditions 

In order to check the effects of the computed mass-loss 
recipe on evolutionary tracks of HB stars, we have com- 
puted two series of models, identical in every way, except 
that one set incorporates the mass-loss recipe, while the 
occurrence of mass loss is neglected in the other. All mod- 
els have been computed using the franec evolutionary 
code (Cassisi & Salaris 1997, Castellani et al. 1997 and 
references therein). As far as the adopted physical input 
parameters as well as the treatment of convection during 
the central He-burning phase are concerned, we refer the 
interested reader to the papers by Cassisi & Salaris (1997) 
and Zoccali et al. (2000). The treatment of outer bound- 
ary conditions is performed, as usual, by adopting a T(t) 
relation (Krishna-Swamy 1966). 

To check the validity of this assumption for hot HB 
stars, we have computed a large grid of model atmospheres 
that provide more accurate descriptions of the thermal 
stratification of the atmospheres of these stars, and inves- 
tigated whether a different treatment of the outer bound- 
ary conditions influences the evolutionary output parame- 
ters. To this end, we computed several HB models for dif- 
ferent metallicities by adopting in one case the T(r) rela- 
tion from Krishna-Swamy, and in the other the boundary 
conditions provided by the more sophisticated model at- 
mosphere computations. The model atmospheres used for 
this test were the non-LTE ISA- WIND models of de Koter 
et al. (1993, 1997) with negligible mass loss, as well as the 
hydrostatic LTE, line-blanketed models of Kurucz (1993). 
In the model atmosphere cases, the connection between 
the atmospheres and the internal structures has been fixed 
at tros = 10. Note that we have verified that the results 
obtained are not affected by this choice of matching point 
between stellar atmosphere and internal regions. Our nu- 
merical experiments clearly showed that the stellar prop- 
erties such as the effective temperature, are not signifi- 
cantly affected by the assumptions made concerning the 

4 http://astro.ic.ac.uk/~jvink/ 
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outer boundary conditions. This is because the T(t) ap- 
proach provides an estimate of the thermal stratification 
in the stellar atmosphere that is in good agreement with 
those provided by more accurate model atmospheres com- 
putations, as can be seen in Fig.||, where we have plotted 
the thermal stratifications provided by the T(t) method, 
and the model atmospheres quoted above, for values of the 
effective temperature and surface gravity suitable for hot 
HB stars. These results imply that canonical theoretical 
predictions of effective temperatures of hot HB stars can 
be considered to be robust. 

5.2. HB Evolution with Mass Loss for metal rich 
clusters 

In order to maximise the effects of mass loss on the HB 
evolution, and given that mass loss increases with stellar 
metallicity, we have computed HB models for a metallicity 
twice solar^. This choice allows a direct comparison with 
the evolutionary computations performed by Yong et al. 
(2000). All HB models have a 1 M Q RGB progenitor with 
an initial chemical composition of Z = 0.04 and Y — 0.34. 
The He core mass of this structure at He ignition is equal 
to 0.466 Mq, while its surface He abundance in the same 
evolutionary phase is of the order of Y = 0.36. The stan- 
dard models, i.e. the ones computed neglecting mass loss, 
have been presented by Bono et al. (1997). 

In the various panels of Fig.^, a comparison between 
standard models and models accounting for mass loss (ac- 
cording to Eq. ||) is presented. When computing the mod- 
els with mass loss, we have accounted for this process along 
the whole evolutionary path, starting from the Zero Age 
Horizontal Branch (ZAHB) until an effective temperature 
of the order of 40 000K (see previous discussion) . This im- 
plies that we are also using our mass loss recipe out of 
its validity range in luminosity. In the case of the most 
massive star, i.e. the coolest one, with a ZAHB location 
below 10 000 K, we are also slightly extrapolating Eq. (||) 
out of its validity range in T c g. 

We would expect that the computed mass-loss rates 
of the order of 10 -12 M yr _1 are too low to alter evolu- 
tionary tracks in a major way. Indeed, as shown in Fig.^J, 
the evolutionary paths of the selected models are not sig- 
nificantly affected by the occurrence of mass loss at the 
computed rates. The less massive, hottest model has lost 
an amount of mass of 2.2 • 10 _4 M Q at the end of the He 
central burning phase, while the coolest model, the one 
with mass equal to 0.530M Q , has lost 2.4 • 10 _4 M Q at the 
end of the same evolutionary phase. The amount of mass 
lost is slightly larger for the more massive model since this 
model is brighter, and the mass-loss rate most strongly 
depends on L. From these evolutionary computations, we 
arrive at the following conclusions: 



5 All evolutionary models computed here as well as addi- 
tional models can be obtained upon request 



— Mass-loss rates due to radiation pressure are not suffi- 
cient as to significantly modify evolutionary tracks for 
HB stars. 

— The mass-loss computations do not provide the high 
rates needed in the Yong et al. (2000) scenario for 
explaining the occurrence of EHB stars in metal-rich 
clusters such as NGC6791. 

— Since the efficiency of radiation-driven wind decreases 
with stellar metallicity, the effects are likely to be even 
smaller for HB stars in more metal-poor clusters. 

5.3. HB Evolution with increasing Mass Loss due to 
radiative levitation 

Before we present the effects of increased mass loss due to 
radiative levitation on HB evolution, we first discuss the 
connection between mass loss and the changes of the sur- 
face chemical abundances due to these physical processes 
of radiative levitation and atomic diffusion. Our calcula- 
tions show that the mass-loss efficiency increases strongly 
when the envelope abundances of heavy elements increase 
as a consequence of radiative levitation. Nevertheless, it 
is also well known that mass loss works as a competing 
process to diffusion by decreasing the efficiency of radia- 
tive levitation in producing large chemical overabundances 

- at least for heavy elements such as silicon. In particu- 
lar, Michaud & Charland (1986) have shown that if mass 
loss increases beyond 10 -14 M Q yr _1 , chemical overabun- 
dances could be wiped out. In more recent times, Unglaub 
& Bues (2001) have investigated the influence of diffusion 
and mass loss on the chemical composition of sdB stars. 
The main outcome of their work was that observed chem- 
ical patterns can only be explained if mass-loss rates are 
in the range 10~ 14 <M(M Q yr _1 ) < 10~ 12 . Higher rates 
would basically prevent the effects of diffusion, whereas 
for lower rates helium would sink in too short time scales 
compared to the typical lifetime of an sdB star. Note that 
our mass-loss predictions fall in the middle of the range of 
the calculations by Unglaub & Bues (2001) for metallici- 
ties typical of GCs with extended blue tails; once radiative 
levitation becomes effective in strongly increasing the stel- 
lar metallicity, mass loss strongly increases, and becomes 
of the same order of magnitude to - or even larger than 

- the upper limit quoted in the work of Unglaub & Bues. 
If these physical processes indeed occur in real stars, they 
should have strong effects on the measured abundances of 
helium and heavy elements in hot HB stars. A detailed 
investigation into this topic is beyond the scope of the 
present work, but we wish to emphasise the further need 
for accurate chemical separation computations that con- 
sistently account for diffusion, radiative levitation, and re- 
liable mass-loss predictions. This, as a function of metal- 
licity and other stellar parameters such as mass, effective 
temperature and luminosity. 

We now turn to the effect of increased mass loss on the 
evolution of HB stars. Given the importance of the metal- 
licity on the line driving efficiency, we maximise the effects 
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Fig. 3. Upper panel: Comparison be- 
tween the atmosphere thermal stratifi- 
cation provided by the T(r) integration 
and by the model atmosphere of Kurucz 
(f 993) for the labelled values of effective 
temperature and gravity. Bottom panel: 
As upper panel, but for the model at- 
mospheres of de Koter et al. (f993, 
f 997), for the labelled values of effective 
temperature and gravity. In all cases, 
a chemical composition of Z = 0.012; 
Y = 0.25 has been adopted. 



of radiative levitation on the surface abundance. This was 
done by simply assuming that starting from the ZAHB 
the surface metallicity is equal to five times or ten times 
the solar value. Note that we do not properly account for 
radiative levitation in the evolutionary code, but assume 
that the diffusion process is a fast and efficient mechanism 
to increase the surface metallicity (Michaud et al. 1985). 
Although this is a rather crude approach, we are only in- 
terested in checking the maximum effect of radiative levi- 
tation on the mass-loss efficiency, and we do therefore not 
account for possible changes in the envelope opacity prop- 
erties due to this metallicity increase. In case we would 
have accounted for these effects, the models would have 
been somewhat fainter and cooler, which would have had 
the effect of slightly decreasing the mass-loss efficiency. In 
other words, our more crude approach has the effect of 
maximising the efficiency of mass loss. 

The additional experiments (also shown in Fig.|J), 
clearly show that even with the assumed increase of the 
surface stellar metallicity, the mass-loss efficiency is still 
too low to affect the HB evolution. Therefore, it is worth 
noting that: 

— It is almost impossible that an in situ mass loss process 
due to radiation pressure alone can affect the HB mass 
distribution in any major way. 

— A relevant parameter in population synthesis studies 
is the maximum mass of HB structures that at the 
end of the central He-burning phase behave as AGB- 
manque stars (Greggio & Renzini 1990). Due to the 



low efficiency of mass loss during the HB phase, this 
critical mass is not affected by the inclusion of our 
mass loss recipe in stellar model computations. 



6. Mass loss, HB rotation rates, & the log g jump 

In this section we discuss the implication of mass loss on 
HB rotational velocities and the "low gravity" problem for 
hot HB stars, which is represented by a jump in log g at 
T off ~11000 K. 

Observational analyses of rotational velocities of HB 
stars show that among the cool group (T e ff < 11 000 K) 
both fast and slow rotators are present (Peterson et al. 
1995), but that for the hot group (T cff > 11 000 K) all stars 
rotate slowly (Behr et al. 2000a, 2000b ; Recio-Blanco et 
al. 2002). We argue that the absence of fast rotators in 
HB stars hotter than 11 000 K can be explained by a stel- 
lar wind set up by radiative levitation of heavy elements 
(Sweigart 2000), which could contribute significantly to 
the removal of angular momentum. Recio-Blanco et al. 
(2002) argue that the mass-loss rate may increase by a 
large factor between 10 000 and 20 000 K due to a change 
in the ionisation state of hot star winds (referring to the 
work by Vink et al. 2000 on massive stars). However, as 
we have shown in Fig. [j], dramatic changes in the mass-loss 
rate over T e g are absent for HB stars. This notwithstand- 
ing, as can be noted from Eq. (||), mass loss increases by 
about 2 dex when the photospheric metal abundance (Z) 
increases by 2 dex. We therefore argue that the striking 
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Fig. 4. The H-R diagram representing the evolution of selected HB models with various initial mass (as labelled) 
whose RGB progenitor mass is equal to 1 M©; with an initial chemical composition of Z = 0.04 ; Y = 0.34. In each 
panel, the short dashed line corresponds to the standard model (neglecting mass loss), while the solid line refers to 
models accounting for mass loss according to Eq. (^). In panels (a) and (c), models accounting for an enhancement of 
the surface metallicity (see text for details) are also plotted (dotted line - Z — 0.1; long-dashed line - Z = 0.2). 



photospheric abundances in hot HB stars, which are most 
likely due to the onset of radiative levitation, may spin 
down the surface velocities of HB stars hotter than 11 000 
K, explaining the absence of fast rotation at these temper- 
atures. Note that firmer and more quantitative conclusions 
can only be achieved by understanding the coupling be- 
tween mass loss and stellar rotation (see Soker & Harpaz 
2000). 

In addition, we question whether the increase in the 
mass- loss rate around 11 000 K invalidates the use of hy- 
drostatic model atmospheres (such as those by Kurucz) 
for hot HB stars. For massive O stars it is a well-known 



fact that neglecting winds in model atmosphere calcula- 
tions causes systematic errors in the derived atmospheric 
parameters, notably log<?. This has even led to a sys- 
tematic discrepancy between masses derived from stellar 
spectra vs. those from evolutionary models, the "mass- 
discrepancy" (Groenewegen & Lamers 1989; Herrero ct 
al. 1992). Indeed, the luminosities for HB stars are much 
lower (having the effect of lowering the mass-loss rates), 
but the stars are also less massive (increasing M), and 
have much smaller radii, which substantially increases the 
more relevant parameter, the mass flux. It is therefore 
not at all obvious that hydrostatic model atmospheres are 
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applicable to these types of objects. Even more so, atmo- 
spheric analyses for hot HB stars have shown that there 
is also a "mass discrepancy" for these objects (Moehler 
et al. 1995), whereas the atmospheric determinations and 
canonical evolutionary models do agree for the cooler HB 
stars. Although the "low gravity" problem for hot HB 
stars may partially be explained by radiative levitation 
of metals (Grundahl et al. 1999; Moehler et al. 2000), it 
has not been solved completely (Moehler 2001), the resid- 
ual discrepancy that is still left between 15 000 < T c s < 
20 000 K is about 0.1 dex in log 5. 

To test whether stellar winds have a noticeable effect 
on stellar spectra of hot HB stars, we compute H7 line 
profiles for these objects using the model atmosphere code 
ISA- wind and the synthetic spectrum code wynspec (de 
Koter et al. 1997). We compare H7 line profiles for a so- 
lar metallicity star with the following stellar parameters: 
T GfT = 17 500 K, M = 0.56M Q , and log(Z/L Q ) = 1.37, cor- 
responding to log g = 4.73. Using Eq. (||), we expect such 
a star to have a mass-loss rate log M (Mqyi^ 1 ) = —11.85, 
and we therefore use this value in the computation of an 
H7 line profile, and compare this line profile to a model 
with negligible mass loss (we choose a rate typical for the 
Sun: logM (Moyr" 1 ) = -14). 

As can be seen in the upper panel of Fig.|| wind 
emission in the model with significant mass loss (dashed 
line) has a noticeable effect on the line wings, which may 
mimic a lower log g from Balmer H7 line profiles. To check 
whether the logy jump could indeed be an artifact of the 
use of hydrostatic model atmospheres, we compare H7 line 
profiles for models with different log g values: log g = 4.73, 
vs. a model with log g = 4.63 (or equivalently a mass of 
M = 0.44Mq), and keep all other parameters fixed. This 
comparison is also shown in Fig.|H| The bottom panel 
shows H7 for the high gravity (solid line) and the low 
gravity (dashed line) model. As expected, Stark effects 
broaden the line wings for the higher gravity model. 

Although the effects in the two panels of Fig.|| are 
small, they yield the same systematic trend. They are so 
similar that it suggests that the neglect of winds in atmo- 
spheric analyses of hot HB stars can mimic too low surface 
gravities. An increase of mass loss due to the increased sur- 
face abundances of hot HB stars may therefore invalidate 
the use of hydrostatic atmospheres for these stars. Note 
that the use of such models is considered robust for the 
cooler HB models, where indeed no discrepancy between 
evolutionary and spectroscopic masses has been reported. 
We argue that the "wind emission" effect - together with 
the metal enhancement due to radiative levitation - is the 
most likely explanation for the observed jump in log g. 

We then arrive at the following scenario for the "zoo" 
of problems in HB morphology around 11000K: due to 
the more stable atmospheres of the hotter HB stars radia- 
tive levitation increases the metal abundances for these 
objects. This subsequently explains: (i) the striking sur- 
face abundances in hot HB stars; (ii) the existence of the 
Stromgren u-jump. It also reduces the logg discrepancy. 
As the abundance spectroscopic determinations indicate 
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Fig. 5. Upper panel: the normalised flux at H7 for a model 
with an expected mass-loss rate of logM(M©yr _1 ) = 
— 11.85 (dashed line), and a model for negligible mass loss, 
i.e. logM(Af Q yr~ 1 )~ —14 (solid line). Bottom panel: a 
comparison for the high gravity (log g — 4.73; solid line) 
and the low gravity (log g = 4.63; dashed line) model. The 
other stellar parameters in both panels are: T e g = 17500 
K, log(L/L©) = 1.37, and solar metallicities. 

that only the hotter HB models are affected by radiative 
levitation, a stellar wind is expected to be set up with 
an increased M by two orders of magnitude. This onset 
of (significant) mass loss can then also help in explaining 
(iii) the absence of fast rotators for HB stars, and (iv) 
the jump in log g. In passing, we finally note that (v) the 
gap observed along the HB of many GCs at (B — V) ~ 
(T c ff ~ 9 000 K) could also be due to an atmospheric phe- 
nomenon related to the chemical peculiarities induced by 
radiative levitation as suggested by Caloi (1999). 

7. Conclusions & Outlook 

In this paper we have, for the first time, computed mass- 
loss rates for HB stars. We have shown that the computed 
rates, as predicted by the most plausible mechanism of ra- 
diation pressure on spectral lines, are too low to produce 
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EHB/sdB stars. This invalidates the scenario outlined by 
Yong et al. (2000) to create these objects by excessive mass 
loss on the HB. We argue, however, that mass loss plays 
a role in the distribution of rotational velocities of hot 
HB stars, and for the so-called logg jump. The mass loss 
recipe derived in this paper is, strictly speaking, only valid 
for HB stars, but as there are hardly any mass-loss predic- 
tions available for low-mass blue stars, the recipe may also 
be applied to: post-HB, AGB-manque, UV-bright stars, 
extreme helium stars, as long as the desired accuracy is 
within a factor of two, and as long as the effective tem- 
peratures are not higher than 40 000 K. 

Although we have proposed a scenario where winds are 
ubiquitous for hot HB stars, and subsequently affect the 
rotational velocities, as well as the atmospheric parame- 
ters (log g), there is still a lot of work to be done. 

First and foremost, spectral evidence for mass loss in 
HB and sdB stars ought to be sought to check whether the 
mass-loss rates, as derived in this paper, indeed occur. 
Diffusion calculations including mass loss for sdB stars 
(Unglaub & Bues 2001) suggest that our derived mass- 
loss rates are reasonable, but this is certainly not a model- 
independent check. 

Second, evolutionary models including rotation (first 
steps have been undertaken by Sills & Pinsonneault 2000) 
and mass loss should be computed to see whether the ab- 
sence of fast rotators for stars hotter than 1 1 000 K can 
indeed be due to the removal of angular momentum due 
to stellar winds. 

Last but not least, systematic atmospheric analyses 
of hot HB stars accounting for the actual surface heavy 
elements distribution and including mass loss should be 
performed to see whether the log g jump is indeed an ar- 
tifact of the adopted hydrostatic model atmospheres. The 
current situation, where evolutionary models are not in 
agreement with the spectral analyses is highly undesir- 
able, as this suggests that current stellar evolution theory 
is not only incapable of producing extreme HB stars, but 
that even "normal" blue HB stars pose a serious problem. 
In other words, a solution to the "low gravity" problem 
for hot HB stars could significantly enhance our current 
understanding of the later phases of stellar evolution. 
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